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This paper proposes groove-like potential structures for the observation of quantum information 
processing by trapped particles. As an illustration the effect of quantum statistics at a 50-50 beam 
splitter is investigated. For non-interacting particles we regain the results known from photon 
experiments, but we have found that particle interactions destroy the perfect bosonic correlations. 
Fermions avoid each other due to the exclusion principle and hence they are far less sensitive to 
particle interactions. For bosons, the behavior can be explained with simple analytic considerations 
which predict a certain amount of universality. This is verified by detailed numerical calculations. 
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The duality between wave and particle aspects is one of the central issues of Quantum Mechanics. Much has been 
0^ ■ made of the particle aspects of photons, but the recent progress in cooling and controlling atomic motion has brought 
forward their wave mechanical behavior in a prominent way. The new field of Atomic Optics has emerged |IJ . 

With modern cooling and trapping techniques, one can envisage controlled motion of atomic particles in structures 
' whose mechanical dimensions match the heterostructures used in electronic circuits. Neutral atoms can be stored in 
^v^j magneto-optical traps, and Hansch and his group has recently shown ]2| that such traps can be made very small, viz 
■ of the order of 10 2 /xm. This requires high precision in the fabrication of the solid structures defining the dimensions of 
t-H | the trap. Modern lithographic technology suggests that such structures could be made even much smaller, and then 
• we can imagine experiments in traps of genuinely microscopic dimensions, where quantum effects would dominate 
the particle dynamics. Hansch has also suggested that such traps could be made into channels and structures, thus 
. providing a tool to design arbitrary devices at the surface of a substrate. 

Similar structures can be constructed by combining charged wires with evanescent wave mirrors ||4|] or magnetic 
Oh, mirrors Such combinations can be used to build up the structures utilized in nano-electronics. The use of wires 
' to guide atomic motion has been investigated by Denschlag and Schmiedmayer Q . Schmiedmayer has also discussed 
the use of such structures to construct quantum dots and quantum wires for atoms . 

Alternative ways to achieve guided motion and possibly controlled interaction between atoms is to utilize hollow 
optical fibers with evanescent waves trapping the atoms to narrow channels at the center of the fiber |^,^| . These can 
eventually be fused to provide couplers similar to those used for optical signal transmission in fibers. Also the pure 
atomic waveguide achievable by the use of hollow laser modes may be used. 

We see these methods as an opening to novel and innovative uses of particle traps. By arranging a network of 
grooves on a surface, we can launch particles (wave packets) into the various inputs of the system, let them propagate 
through the device and interact with its structures and each other. This may well provide an opportunity to design 
quantum apparatuses, process information and perform computations. The advantage is that both the structures 
and the input states are easy to control in an atomic environment. An equivalent point of view is expressed by 
Schmiedmayer in Ref. @. 

A next step in the experimental progress would be to observe the quantum character of atoms (or possibly ions). 
An essential quantum characteristic of particles is their statistical behavior. The difference between bosons and 
fermions manifests itself dramatically in many situations. Optical networks can be fed by a few photons only, and 
their quantum aspects have been utilized in experiments ranging from secure communication to tests of fundamental 
issues. Recently Zeilinger and his group have tested the behavior of two-photon states at beam splitters. Using 
the overall symmetry properties of the states, they have been able to display both symmetric and antisymmetric 
behavior. 

Similar experiments are in principle possible with electrons. In nanostructures, one can fabricate the devices 
simulating optical components, but it is far less trivial to launch single conduction electrons in well controlled states. 
Yamamoto's group, however, has been able to show quantum correlations in an experiment which is the analogue of 
a beam splitter for photons |y]]. 

In this paper we give an example of the multiparticle effects observable when particle states are launched along 
potential grooves on a surface. The specific phenomenon singled out for investigation is the effect of particle statistics 
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at a beam-splitter-like coupling device. The corresponding effect with photons is described in Sec. II as a motivation. 
In Sec. Ill we present the details of the model chosen and a simplified analytic treatment demonstrating the main 
features expected of this model. In Sec. IV we carry through a numerical analysis of the situation, for one particle 
as a two-dimensional propagation problem, but for two particles in a paraxial approximation. For non-interacting 
particles, the expected behavior is found, but when particle interactions are added, the boson behavior is changed. 
For fermions the exclusion principle make them essentially insensitive to the interaction. An unexpected feature is 
found: the sign of the interaction is irrelevant for the effect. In Sec. V this is explained within our simple analytic 
model, and, as a consequence, a certain universality is proposed: When the interaction strength over the tunneling 
frequency becomes of the order of V3, the noninteracting bosonic behavior is essentialy destroyed. This is verified 
by numerical calculations, reported in Fig. 13. Finally Sec. VI presents a discussion of parameter ranges in real 
materials, where our effects may be observable, and summarizes our conclusions. 



II. MOTIVATION 



In order to show the opportunities offered by atomic networks, we investigate the manifestations of quantum 
statistics on an experiment emulating the behavior of photons in beam splitters. This is a straightforward approach, 
which enables us to display the potentialities and limitations of such treatments. 

Our work has been motivated by the statistics displayed by a 50-50 beam splitter, which has been used in the 
experiments by the Zeilinger group ]l0| . When two particles are directed into the beam splitter in the incoming 
modes in Fig. [lL they are piloted into the outgoing modes according to the beam splitter relations 
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see Ref. |0|. When one particle is directed into each incoming channel, the state is 

l*) = 4XJo>, 



(2) 



where |0) is the vacuum state. Without assuming anything about the statistics of the incoming particles, we can 
express the state (0) in terms of the outgoing states by inverting the relation (pi) as 



(«Lt) 2 + (&L)' 



1 h> 

out ! out 



|0>. 



(3) 



From this follows that boson statistics gives 



I*) = (\n a ,out = 2, n 6 , out = 0) + |?Vout = 0, n Mut = 2)) 
the particles emerge together at either output. For fermions we have 

l*> = oL&L|o> 

= |?Vout = !,™&,out = 1), 



(4) 



(5) 



and they always remain separated. 

Weihs et al. [1(J have been able to verify these properties experimentally using photons. As the requirements of 
quantum statistics refer only to the total wave functions, they have been able to realize both the symmetric and the 
antisymmetric case, thus offering the behavior of both bosons and fermions. 

Photons are ideal for experiments, they do not interact mutually and they propagate essentially undisturbed in 
vacuum. As models for quantum systems, they have the drawback that they cannot be localized, their wave packets 
are of rather elusive character, and the influence of particle interactions cannot be established. Thus we have chosen 
to discuss the propagation of massive particles through beam-splitter-like structures. As explained above, such 
experiments may be performed with atoms or electrons in traps of microscopic dimensions. We can thus investigate 
the propagation of wave packets through these structures, explore the role of quantum statistics and switch on and 
off the particle interaction at will. 
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III. THE MODEL 



We consider particles moving in potential wells which form grooves over a two-dimensional surface. These may 
cross or couple by tunneling when approaching each other, thus forming a network of potential channels emulating a 
linear optical system. 

Here we consider two separate channels which run parallel for z — > ±00 and approach each other in the a;— direction, 
as shown in Fig. |2[ For simplicity we construct the potential from two harmonic oscillators 

U ± (x)= X -mu?{x± X -d) 2 , (6) 

so that a double well potential can be obtained by writing 

z) = U+(x,z)U-(x,z) 
' U+(x, z) + U-(x, z) 

1 2 (x + ^d(z)) 2 (x- \d{z)) 2 

— — JJIU) — 

2 ' (x+\d{z)) 2 + {x~\d{z)) 2 ' 

If we now choose d{z) in a suitable manner, we can achieve the potential behavior shown in Fig. |^. Note that at 
the minima, the potential U(x, z) essentially follows the shape of the smaller potential U±. 

We consider a wave packet sitting stationary near the bottom of one well at z = 0, where the distance between the 
wells is at its minimum do . The particle can then tunnel across the barrier with the rate 



T ~ exp 



y / 2mU(x 1 0)dx 



exp 



K y/U(0,0)do] , (8) 



where K is some constant. From Eq. (Q) we see that £7(0,0) oc dq so that we expect 

logT n'dl + const. (9) 

In order to acquire a heuristic understanding of the physics involved in the coupling of the grooves at z = 0, we 
look at the lowest eigenfunctions of the double well potential. These are expected to be symmetric, ips, with energy 
Es, and antisymmetric, ipA, with energy Ea, as shown in Fig. ||. We have Ea > E$ and hence we write 

e a = e + nn 

E s =E - hfl, (10) 

where 20 is the tunneling frequency. 

Using the eigenstates we form the localized states 

<Pl = ^= (i>s + i/>a) 

= -J= (ipS ~ i>A) » (11) 

where the subscripts L (R) denote left (right) localization. 

We can easily integrate the time evolution by using the energy eigenstates. If we now assume that we start from 
Pi at time t — 0, then 

*(*) = exp(-iHt/h) ip L 

= 4= exp (-iEt/h) (e mt ^ s + e~ mt ^ A ) (12) 
= exp (— iEt/h) (cos fit ip l + i sin fit ipn) ■ 
This displays the expected flipping back and forth between the two wells. For 

Clt = - A (13) 
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the coupling performs the action of a 50-50 beam splitter. 

We now move to consider the action of such a potential configuration on a two particle initial state. We first choose 
the bosonic one 

^ = i= (v>l(1)M2) + ^(2)M1)) , (14) 
where the argument denotes the coordinates of the particle. This can be expressed as 

* B = -L (^(1)^(2) - ^a(1)Vu(2)) , (15) 
which can be evolved in time straightforwardly to give 

exp (-iHt /h) *o=^l cxp (-i2Et /h) (e i20t ° Vs(l)^s(2) - e^ 2nt ^ A (l)^ A {2)) 

= -±= exp (-iiEto/h) (<p L (l)<p L (2) + <Pr(2)<pr(1)) . (16) 

As we see, the bosonic two particle state works as in Eq. (||): both particles emerge together. 
In the fermionic case we have 

\£f = (^(1)^(2) - ^(2)^(1)) 

= ^=(^a(1)Vs(2)-Vs(1)V'a(2)). (17) 

Because both states ip A ips and V'sV'A evolve with the energy 2E, remains uncoupled to other states. Thus the 
fcrmions emerge at separate exit channels as expected. 

IV. NUMERICAL WORK 
A. The Schrodinger equation 

The Schrodinger equation in the two-dimensional system is of the form 

* 9 T / \ r h 2 ( d 2 d 2 
th—*(x,z,t)= I - — 



2m V dx 2 ^~ dz 2 



V(x,z,t). (18) 



As a preparation for the numerical work, we introduce the scaling parameters r and £ giving the dimensionless 
variables 

x = x/£ 

z = z/i (19) 
t = t/r 

p = Tp/mt;. 

We apply this to the one-dimensional oscillator Hamiltonian 

H=£- + \mu?x 2 (20) 
2m 2 v ' 

and find the Schrodinger equation 

if , ~- , , 

mejdt \2 2 

The dimensionless oscillator frequency is given by uj — lot . This shows that choosing the scaling units suitably, we 
can tune the effective dimensionless Planck constant 



i(^)U=(£ + &#)*- (2D 
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To check the consistency of this we calculate 



&a = (l) (^) 

This gives us a way of controlling the quantum effects in the numerical calculations. 
In our numerical calculations we employ the split operator method Jl3| 

exp [-i(T + U)At/h] re exp [-iT At/ft] exp [-iUAt/h] . 

The corrections to this are given by 



At 2 



At C 2 ' 



ij) + pa; 



(22) 



(23) 



(24) 



(25) 



In order to achieve satisfactory numerical accuracy, this should not be too large; in our calculations, with At = 0.001, 
ui = 30 and h = 6 the expectation value of expression (^5|) is of the order of 10~ 4 . Decreasing At or the grid spacing 
has been found not to change our results significantly. 

In the following discussion, we use the scaled variables, but for simplicity, we do not indicate this in the notation. 
Whenever variables are assigned dimensionless values, these refer to the scaled versions. 

In order to achieve beam splitter operation, we let the distance between the potential wells vary in the following 
way 



d(z) = 2 + d 



cosh(z/?y) ' 



(26) 



which inserted into Eq. ([?]) gives a potential surface as shown in Fig. ||. To test its operation as beam splitter, we let 
a wave packet approach the coupling region in one of the channels, and follow its progress through the intersection 
numerically as a two-dimensional problem. The result is shown in Fig. |]. We see that the parameters chosen lead 
to ideal 50-50 splitting of the incoming wave packet. The progress of the wave packet through the interaction region 
is steady and nearly uniform, and no backscattering is observed. This suggest simplifying the situation so that the 
motion in the z-direction is replaced by a constant velocity, and the full quantum problem is computed only in the 
x-direction. If the wave packet is long enough in the z-direction, its velocity is well defined, and this should be a good 
approximation. 

The implementation of such a paraxial approach becomes imperative when we want to put two particles into the 
structure. The full two-dimensional integration would require the treatment of four degrees of freedom, which is 
demanding on the computer resources. With the paraxial approximation, two particles can be treated by a two- 
dimensional numerical approach, which is within the resources available. 

To introduce the paraxial approximation, we perform a Galilean transformation of the wave function to a co- moving 
frame 



ip{x, z, t) = <p(x, t) exp 



Poz 



P& 
2m 



where 



m 



(27) 



(28) 



The initial momentum in the z-direction is denoted by po- The new wave function is found to obey the Schrodinger 
equation 



,- 9 , 



V- 


( d 2 


d 2 \ 


2m 


\dx 2 ~ 


h 0?V 



m 



(29) 



For a well defined momentum po, the wave packet is very broad in the ^-direction and its derivatives with respect to 
<; may be neglected. The corresponding degree of freedom disappears, and it is replaced by a potential sweeping by 
with velocity po/m. This is what we call the paraxial approximation. 
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Taking the parameters from the integration in Fig. ^, we can obtain the beam splitting operation also in the 
paraxial approximation as shown in Fig. a. The transfer of the wave packet from one well to the linear superposition 
is shown in Fig. [fj. This proves that the potential configuration works exactly as in the analytic result (|lj) . 

Now the integration is one-dimensional for a single particle, and it is easy to investigate the tunneling probability 
as a function of the parameters. In Fig. ^ we display the transition rate T as a function of the parameter cIq, which 
controls the coupling between the wells. For small values, do < 3, we are in a coherent flipping region; the wave 
packet is transferred back and forth between the wells and resonant transmission occurs. For larger values, d§ > 3, 
the analytic estimate of a logaritmic dependence in Eq. (^|) is seen to hold approximately. Our calculations work at 
d = 1.8903, which gives T = 1/2. 



B. Effects of quantum statistics 

We can now integrate the propagation of a two-particle wave function by choosing the initial state to be combinations 



of 



-— (x±(l + -d ) 



(30) 



where the + (— ) refers to the particle entering in the left (right) channel. For bosons, this is used in the combination 
(14) and integrated in the potential ([?]), where the z— dependence is replaced by a t— dependence according to Eq. 
(29). The result is shown in Fig. ||. At t = —10, the bosons enter symmetrically in the two input channels, i.e. they 
have different signs for their coordinates. After being mixed at time t — 0, they emerge together with equal strength 
at both output channels, i.e. their coordinates have the same sign. This result fully reproduces the behavior expected 
from bosons at a 50-50 beam splitter. 

We can, however, also test the fermionic case by using the state (jl?]) as the initial one. The result is shown in Fig. 
|9|. Near t — 0, the wave packets follow the potential wells, but they remain separated and emerge at different outputs 
as expected. Fermions do not like to travel together. 

We have thus been able to verify the properties of a 50-50 beam splitter on massive particles represented by wave 
packets travelling in potential structures. The calculations in Figs. (|^) and (^|) do not, however, include any particle 
interactions. We can now proceed to include these, and evaluate their effect on the manifestations of quantum 
statistics. 

When we introduce the interaction, we have to decide which type of physical system we have in mind. Conduction 
electrons or ions interact through the Coulomb force whereas neutral atom interactions may be described by a force 
of the Lcnnard- Jones form. In both cases, the interaction is singular at the origin, and it has to be regularized there. 
We do this by introducing the variable 



(31) 



This makes the interaction energy finite when the particles overlap, but does not affect the main part of our argument 
in other ways. 

With this notation the Coulomb interaction is written 



V c (r) 



V 



(32) 



and the Lcnnard- Jones interaction is 



V LJ (r) = V Q 



b\ 12 fb^ 



(33) 



this contains the additional range parameter b. In both cases, the strength of the interaction is regulated by Vo- 
lt is straightforward to integrate the Schrodinger equation with the two-particle interaction included, and look how 
its increase affects the correlations imposed by quantum statistics. For fermions, the effect is essentially not seen in 
the parameter ranges we are able to cover. As seen from Fig. ^, the particles never really approach each other, and 
they remain separated due to their quantum statistics for all times; the interaction does not affect them. 

For bosons the effect is different. We have investigated their behavior for a range of interaction parameters and 
find that an increase in the interaction does destroy the simple behaviors found for noninteracting ones. The result of 



G 



such an integration is found in Fig. 10. Compared with Fig. |8J this shows that now the particles appear in separate 
channels with about the same probability as in the same channels. Thus the statistical effect has been destroyed. Fig. 
[TT| shows how the result emerges during the time evolution; the system tries to achieve the ideal case, but settles to 
the final state observed in Fig. [To[ 

Figure ^ shows how an increase in the interaction strength destroys the ideal behavior. This is drawn using a 
Lcnnard-Jones potential, but the behavior is similar for other cases we have investigated. 

One unexpected feature emerged from our calculations: The destruction of the ideal behavior turned out to be 
independent of the sign of the interaction. One may have expected an attractive interaction between the bosons to 
enhance their tendency to appear together, but this turned out not to be the case. In order to understand this feature 



we have to turn to our simple analytic argument in Sec. HI and investigate the interplay between two-particle states 
and the interaction. 

V. STATISTICS VERSUS INTERACTIONS 

In order to introduce the particle interaction, we choose a convenient basis for the two-particle states. As the first 
component we choose the bosonic state ( |l4| ) 

«i = ^= (^(1)^(2) + V l{2)<pr{1)) . (34) 

In addition, there are two more bosonic states, where the particles enter in the same channels 

u 2 = -= {<Pl(1)<Pl(2) + <pn(2)<p R {l)) (35) 



V2 
1 

V2 



-4 (^(1)^(2) +^(1)^(2)) (36) 



and 



«3 = -j= (<PlO)M*) - M2)¥*(l)) (37) 
= -^(V>s(l)M2) + Ml)V>s(2)). 



As the last component we have the fermionic basis function ( |17| ) 

u± = -L {<p L (l)<p R (2) - ¥> £ (2)p fl (l)) . (38) 

Together the functions {ui} form a complete Bell state basis for the problem. They are also convenient for the 
introduction of particle interactions. In the states u± and u 4 the wave functions overlap only little, and the effect of 
the interaction is small. For the states u 2 and W3 they sit on top of each other and feel the interaction strongly. We 
treat this in a Hubbard-like fashion by saying that the energy of the latter states is changed by the value 2V oc Vq . 
Because the overlap between ip^ and ipu is small, we have {u\, Vu\) « (w 4 , Vua) and (u 2 , Vu 2 ) ~ (113, VU3), and we 
can use the definition 

2V = X - [(ua, Vu 2 ) + (u 3 , Vu 3 ) - («i, Vm) - (u 4 , VU4)] (39) 



ip(x) 2 (p(y) 2 V(\x - y\)dxdy - I I tp(x) 2 (p(y - d ) 2 V(\x - y\)dxdy, 



where V is either a Coulomb or Lennard- Jones interaction. The first terms give the effective interaction energy when 
both particles sit in the same potential groove. In the states u\ and U4 both grooves are occupied and because they 
are at their closest at z = 0, we subtract the mutual interaction energy across the separating barrier to obtain the 
pure local interaction energy. The second line results if we approximate both ipi, and tpn by a Gaussian if with the 
same width. From Eq. ( Jig ) we see that only the states U\ and u 2 are coupled by the tunneling rate 2f2. Thus if wc 
express the state by 

* = a\U\ + a 2 u 2 + 03^3 + a 4 ii 4 , (40) 
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the state vector [01,02,03,04] evolves with the Hamiltonian 



2E -2hQ 

-2HQ 2E + 2V 

2E + 2V 

2E 



= 2E + V- 



-V 


-2hn 








-2Mi 


V 














V 














-V 



(41) 



The constant part does not affect the coupling between the states, and the amplitudes 03 and 04 decouple. The 
remaining ones flip at the effective rate 



n cS = iyW + v 2 /n 2 . 



(42) 



This result shows that the new parameter replacing fl in Eq. ( |io| ) is fl c s, implying that the perfect boson behavior 
is expected to be destroyed for 



V 

2htt 



V3 = 1.73. 



(43) 



This is in approximate agreement with the result shown in Fig. [lj. By inspecting Fig. O , we can also verify that 
the flipping does occur faster when we switch on the interaction, as expected from Eq. ([42^ 

In the simplified analytic treatment, the only influence of the potential was through its strength \V\. Hence we 
expect the results to scale with the parameter (|V r |/2?ii7) where 2M7 = Ea — E$- The probability to emerge in the 
same output channels should essentially depend on this only; a certain universality is expected. 

In Fig. [l3|, we have plotted this probability for a variety of potentials including both Coulomb and Lennard- Jones 
ones. As we can see, the behavior is very similar, at (|V|/2?Lf2) w 1.7 the probability has decreased to less than 10% 
in agreement with our expectation. This verifies the degree of universality achieved. For comparison, we also used the 
simple analytic theory to obtain the points along the curve. In this treatment, 2fl was assumed to be constant during 
some finite coupling time t, according to Eq. (|l3|). By inspecting Fig. ||, we conclude that t should be of the order 
of unity. Here 2hQ was chosen to be 8, and the time evolution in the subspace {1*1,1(2} was calculated. This agrees 
best with the numerical results for small V] for larger values of V the simple analytic treatment becomes invalid. 



VI. DISCUSSION AND CONCLUSIONS 



The actual values of £, r, m and w depend on the physical system we have in mind. Our calculations are carried 
out at ui = 30 and h = 6; the momentum in the z-direction of the incoming wave packet is p z — 30 or p z = 1000. If 
we consider an atomic beam splitter for Rb atoms, setting the length scale £ to 100 nm corresponds to a time scale 
t of 80 /is according to Eq. (^|]). A displacement of 400 nm from the center of one valley in the x-direction gives 
a potential energy increase of roughly 10 mK, i.e. a transverse velocity of 0.15 m/s. This is to be compared with a 
typical height of the confining potential in a hollow optical fiber, a few tens of mK ^|9). Taking p z = 1000 yields a 
beam velocity of 1.3 m/s in the z-direction. 

If we consider a mesoscopic electron beam splitter built on GaAs, £ could be of the order of 40 nm, which means 
that the closest distance do between the valleys is 80 nm. This would correspond to r — 6 ps, i.e. the electron goes 
through the device in a few picoseconds. With p z — 30, the kinetic energy of the electron due to the motion in the 
z-direction would be of the order of 0.01 eV. A displacement of 100 nm from the center of one valley in the ^-direction 
corresponds to a potential energy increase of roughly 0.05 eV, in comparison with the bandgap in GaAs, 0.115 eV. 

The parameter ranges chosen in our illustrative computations may not be experimentally optimal, but they indicate 
that the effects are not totally outside the range of real systems. Even if our calculations are based on a rather simplified 
model, we find them suggesting effects possible in realistic setups. The main problem is to prepare the appropriate 
quantum states, launch them into the structures and retain their quantum coherence during the interaction. With 
atomic cooling and trapping techniques, this may be feasible in the near future. For electrons the possibility to retain 
quantum coherence is still an open question. 

We have chosen to discuss the straightforward question of particle statistics at a beam splitter. This is a simple 
situation, which, however, presents genuine quantum features. For information processing and quantum logic slightly 
more complicated networks are needed. Simple gate operation can be achieved along the lines described in Ref. [jl4j , 
which was formulated in terms of conduction electrons, but similar situations can be envisaged for interacting atoms. 

In this paper we suggest an analytic model, which can be used to analyze the behavior of particle networks confined 
to potential grooves in two-dimensional structures as discussed also by Schmiedmayer 0. For complicated situations, 
the full numerical treatment rapidly becomes intractable and simplified qualitative tools are needed. We hope to have 
contributed to the development of such methods in the present paper. 
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FIG. 1. Schematic drawing of a beam splitter. The incoming modes ai n ,6i n are piloted into the outgoing modes a ou t,&out 
according to the beam splitter relations. 
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FIG. 2. A contour plot of the beam splitter potential U{x,lz). The two channels approach each other at 



variables x and z used in the numerical calculations are related to x and z according to Eq 
oscillator frequency is u = 30. The distance d(z) between the two valleys is chosen as in Eq. (|2(i ) 
Two cross sections of the potential at x — and x = 2.5 are also shown. 
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FIG. 3. Symmetric and antisymmetric eigenfunctions ips and ipA of the double well. 
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FIG. 4. Time evolution of a two-dimensional Gaussian wave packet propagating through the 50-50 beam splitter. The frame 
shown is centered around the wave packet moving in the z-direction. The dashed lines represent the minima of the potential. 
The potential is given by the scaled version of Eq. (Q) with d(z) as in Eq. (|2(j) where r\ = 30. The scaled oscillator frequency 
is lu — 30; h — 6. 
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FIG. 5. Time evolution of a wave packet propagating through the 50-50 beam splitter in the paraxial approximation. The 
wave packet is incident in one channel and splits into two equal parts. Parameters are as in Fig. ^, the only difference is 
that in the potential U(x,z), z is replaced by tpo/m according to the paraxial approximation. The initial wave packet is 
<p°{x) = N exp [-§- (x - (1 + do/2)) 2 ] . 
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FIG. 6. The probabilities, as a function of time, to find the particle of Fig. |B| in the right (full line) and left (dashed line) 
valley of the beam splitter. The particle is incident at the right input and emerges with equal probability at the left and right 
outputs. Parameters as in Fig. U. 
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FIG. 7. Tunneling probability as a function of the square of the minimum distance between the valleys, do- The point of 
50-50 beam splitter operation, do = 1.8903, is indicated with dashed lines. For small values of do, the wave packet is transferred 
back and forth between the wells and resonant transmission occurs. For larger values, do > 3, the relation Q is seen to hold 
approximately. 
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FIG. 8. Two bosons propagating through the beam splitter. Snapshots of the two-particle wave-function |^| 2 at different 
times are shown. On the horizontal axes we see the coordinate of particle 1; the vertical axes refer to particle 2. The two 
particles are seen incident in different input channels; they mix around t = and finally exit together. Parameters as in Fig. n. 
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FIG. 9. Two fermions propagating through the beam splitter. Snapshots of the two-particle wave function \ty\ 2 at different 
times are shown. On the horizontal axes we see the coordinate of particle 1; the vertical axes refer to particle 2. The two 
particles enter and exit in different input and output channels. Parameters as in Fig. |f| 
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FIG. 10. Two bosons propagating through the beam splitter; interaction is included. Snapshots of the two-particle wave 
function \^\ 2 at different times are shown. On the horizontal axes we see the coordinate of particle 1; the vertical axes refer 
to particle 2. The two particles enter in different input channels and mix around t = 0. With no interaction the two bosons 
were always emerging together, as in Fig. H. Now there is a finite probability also for the bosons to emerge in separate output 
channels. The interaction is given by Eq. (pi|), with Vq = 50 and e = 1. Other parameters as in Fig. |E| 
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FIG. 11. The probabilities, as functions of time, to find the two bosons of Fig. ^ in different valleys (full line) and in the 
same valley (dashed line). 
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FIG. 12. The probabilities, as functions of interaction strength, for two bosons to emerge at different (full line) and same 
(dashed line) output channels. The Lennard- Jones interaction is given by Eq. (^3|), with b — 0.25 and e = 0.2. Parameters for 
the potential are as in Fig. H 
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FIG. 13. The probability, as a function of interaction strength, for two bosons to emerge at the same output channel for 
several different types of interaction (full lines). We see that the perfect bosonic behavior is destroyed for V/2h£l ~ \/3 = 1.73. 
The parameter ranges are e — 0.1 to 1 for the Coulomb interaction and 6 = 0.25 to 0.5, e = 0.2 to 0.35 for the Lennard- Jones 
interaction. The numerical results are compared with that of the analytic treatment (diamonds). In this model 2hQ was chosen 
to be 8, and the time evolution in the subspace {u\,U2} was calculated. 
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